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Abstract — In this paper, a parametric cliannel estimation 
algoritiim applicable to uplinks of orthogonal frequency division 
multiple access (OFDMA) systems whose subcarriers are pseudo- 
randomly allocated is proposed. By exploiting pilot hopping, 
estimation of signal parameters via rotational invariance tech- 
nique (ESPRIT) is employed to estimate the path delays of the 
sparse multipath fading channel. From the delay information, a 
channel interpolator utilizing global pilots, which can be irregular 
distributed, is derived to estimate the channel state information 
on the desired tones. Moreover, a simple method of estimating the 
time correlation of the channel taps is introduced and integrated 
in the proposed algorithm. Simulation results demonstrate that 
the proposed algorithm outperforms the local linear channel 
interpolator within a wide range of SNR's and Doppler's. 

Index Terms — Parametric channel estimation. Uplink OFDMA, 
Pilot hopping, ESPRIT. 

L Introduction 

Recently, the pseudo-random resource allocation scheme is 
adopted by wireless standards such as IEEE 802. 16e wireless 
MAN standard [1] to provide frequency diversity and co- 
channel interference averaging in orthogonal frequency divi- 
sion multiple access (OFDMA) systems. However, the price 
paid for the flexibility of such complex allocation scheme 
is the increased difficulty in estimating the channel impulse 
response (CIR) or channel frequency response (CFR) of the 
uplink. Since the pilot tones allocated to a certain user 
are pseudo-randomly interleaved with the others, traditional 
channel interpolators [2] [3] [4] using global pilots are im- 
practicable, while the local linear interpolators suffer from 
significant performance degradation at high SNR regime where 
the estimation error floor shows itself. 

In order to overcome the irregular pilot pattern, [5] has 
proposed a parametric channel estimation [6] scheme which 
can greatly reduce the channel estimation error for sparse 
multipath channels. [5] has applied estimation of signal param- 
eters via rotational invariance technique (ESPRIT) [7] which 
exploits the shift-invariance property of uplink tile structure 
of IEEE 802. 16e to estimate and track multipath delay spread 
of the received signal. With the multipath delays, the channel 
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frequency correlation information can be derived to enable the 
sinc-function interpolator to estimate the CFR on the desired 
tones. 

The downside of [5] is its dependence on the special tile 
structure which introduces the shift-invariance property — 
the base of ESPRIT, however this symmetric structure may 
be corrupted. For example, when uplink virtual multiple- 
input multiple-output (MIMO) is active, pilots in each tile are 
divided into two exclusive subsets and allocated to two users, 
respectively. For any one, the pilots in the subset lose the shift- 
invariance property. As a result, [5] fails under uplink virtual 
MIMO. 

In this paper, the strict restriction of the special symmetry 
structure of pilot pattern is relaxed, in stead, pilot hopping is 
exploited to enable the parametric channel estimation under 
arbitrarily irregular pilot distribution. By combining two vec- 
tors of pilots of two contiguous OFDMA symbols, the shift- 
invariance property is recovered. The sample auto-correlation 
matrix of the combined pilot vectors is formed, and from 
which, time correlation of the multipath fading channel is 
estimated to compensate the time variance induced by Doppler 
Minimum description length (MDL) [8] is adopted to identify 
the number of significant paths. With the delay information, a 
global channel interpolator is introduced to estimate the CFR 
on the desired tones. 

The rest of this paper is organized as follows. Section II 
introduces the channel model and the OFDMA system model. 
The proposed channel estimation algorithm is presented in 
Section III. Simulation results and analysis are provided in 
Section IV. Finally, Section V concludes the paper. 

Notation: Boldface letters denote matrices or column vec- 
tors. tr{-), (•)*, (■)^, (■)^, and || • ||f denote trace, conjugate, 
conjugate transposition, Moore-Penrose pseudo inverse, and 
Frobenius norm, respectively. E{-) represents the expectation 
of a stochastic process. denotes the (i,j)-th element of a 
matrix. span{A) denotes the column space spanned by A. 

II. System Model 

Consider the uplink of an OFDMA system with a bandwidth 
of BW = l/T Hz (T is the sampHng period). We use N to 
denote the total number of tones and Nused the total number 
of useful tones, including data and pilot tones. A cyclic prefix 



(CP) of length Lcp is inserted before each OFDMA symbol 
to eliminate inter-block interference. Thus the whole symbol 
duration is = (TV + Lcp)T. 

Now we consider a certain user is scheduled to transmit over 
M (M < Nused) tones per symbol. The indexes of these tones 
are collected in a set, denoted as Io{n), which can be changed 
along the time. Regardless of CP, the transmitted signal of this 
user can be expressed in matrix form as x(n) = F^s(n), 
where s(n) G qNxi Jj^^ non-zero elements at Xo(n), and 
F e (jNxN balanced Fourier transform matrix with 

the (m, n)-th entry ^e-J^^rnn/N^ 

The discrete complex baseband representation of a mobile 
wireless CIR of length L can be described by [9] 



L-l 



h{'n',T) = ^ Ji{n)6 (t - n) 



(1) 
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where r; e 7?. is the non-sample-spaced delay of the l- 
th path normalized by the sampling period T, and 7((n) 
is the corresponding complex amplitude. Due to the motion 
of mobile stations (MS's), 7i(n)'s are wide-sense stationary 
(WSS) narrowband complex Gaussian processes, which are 
uncorrelated among different paths based on the assumption 
of uncorrelated scattering (US). Furthermore we assume delays 
of multipath are static during the process of estimation, as they 
are varying very slowly compared with the complex amplitude 
of paths [10]. We collect t;'s into a set, denoted as X^. 

In addition. We assume that 7((r7,)'s have the same normal- 
ized correlation function rt{m) for all /'s, hence 



(2) 



where af is the power of the l-th path. 

Assuming a sufficient CP, i.e., Lcp > tl-i, the received 
signal of the base station (BS) for the given user, denoted as 
y(n) S C^^^, can be described as 



y(n) = X(n)H(n) + v(n) 



where X(n) G qNxn ^ diagonal matrix with non-zero 
elements drawn from x(n), H(n) e C^^^ is CFR experienced 
by the n-th OFDMA symbol, and is expressed as 



H(n) = F^h(n) 



(4) 



where h(n) = [h{n,To), ...,h{n,TL-i)f , and F^ e C^^-^ 
is the non-balanced Fourier transform matrix with the (fc, l)- 
th entry g-J^-^kn/N^ ^^^^^ ^ qNxi zero-mean 

complex Gaussian noise vector with variance ct^, i.e., v(n) ~ 

III. Channel Estimation 

Denote the number of pilots for the given user as _P < 
M, the indexes of pilots are collected into a set denoted as 
2p(n) C loin). At the BS, pilots are firstly extracted after 
Fourier transformation to perform least-squared (LS) channel 
estimation in frequency domain, i.e., 

Hp{n) = X-\n)yp{n) = Hp{n) + Wp{n) (5) 



where, Xp(n) G C is the pre-known pilot diagonal matrix, 
Hp(n) e C^'^^ is the estimated CFR, and Wp(n) G C^^'^ is 
the noise vector expressed as 

Wp(n) = Xp ^(n)vp(n) (6) 

and Wp(n) ~ CA/'(0, a^Ip) when X^X^ = Ip for PSK 
modulated pilots. 

In the following section, we will introduce a simple pilot 
hopping scheme based on which the parametric channel esti- 
mation can be carried out. 

A. Pilot Hopping Pattern 

Pilot hopping is made in two phases: the irmer hopping 
between two contiguous OFDMA symbols and the outer 
hopping superposed over the inner one. 

First, the inner hopping is carried out, which shifts all pilots 
of odd OFDMA symbols by a pre-defined offset denoted 



as V. Therefore, after that, X 



V, where X, 



and Xp are indexes of pilots of odd and even OFDMA 
symbols, respectively. Second, the outer hopping is performed 
by shifting the pilots of both even and odd symbols with the 
same offset. 

For the proposed channel estimation algorithm, the inner 
hopping is indispensable, while the outer hopping is optional 
— although it can accelerate the estimation of the auto- 
correlation matrix of CFR by simulating more channel fading 
at the receiver through including phase rotation [11]. For 
simplicity, the outer hopping is not applied in this paper. Later 
on, the shift-invariance structure is exploited from the inner 
pilot hopping scheme. 

B. Estimation of Multipath Delays 

Stack two contiguous estimated CFR vectors into one with 
the even on the odd, i.e.. 



(3) Hco„,p(n) 



Hp(2n) 
Hp(2n + 1) 



+ Weo„,p(n) 



(7) 

where 'Rcon,p{n),'H.con,p{n) G C'^^^^ are the concatenated 
CFR vector and its estimation, respectively, and ^con,p{n) G 
q2Pxi jg jj^g corresponding noise vector, and Wg 

From (4), we have 



^ con^p 



Hp(2n) = Fp,eX2n) 
Hp(2n + 1)= Fp,odh(2n + 1) 



(8) 
(9) 



where Fp^e„,Fp^orf € C^^^ are submatrices drawn from F^ 
with rows indexes belonged to Tp^ev and Tp^od, respectively. 
P > L is required so that Fp^^y and Fp^od are of full column 
rank. As Ip^od — Ip,ev = it is straightforward that 



■ p^od 



F * 



(10) 



where * G C^^^ is a diagonal phase- twisted matrix with 

the l-th diagonal element = e^^^^"^'^^. Apparently, * 
contains the multipath delay information rfs as expected. 



According to (8)-(10), (7) can be rewritten into 



H 



con^p (^) 

" F 







F * 



h(2n) 
h(2n+ 1) 



+ Wc 



(11) 



Consequently, the auto-correlation matrix of the concate- 
nated estimation of CFR is given in (12) shown at the 

bottom of the next page, where Jli, {rri) G C^^^ denotes 
the auto-correlation matrix of the CIR vector with the lag m. 
From the WSSUS assumption and (2), the (i,j)-th element 
[R/,(m)].^- = afrt{m)S{i — j), where S{-) is the Kronecker 
function. Therefore, R?i(m) can be written as 



Hhim) = rtim)Kh{0) 



(13) 



where R?i(0) is a diagonal matrix, and the l-ih diagonal 
element [R;i(0)]; ; = af. 

Discarding the noise component temporally, the useful com- 
ponent on the right-hand side of (12) is rewritten into 



R, 



■con,p 



'^p,ev^h{^)'^ p^ev 

r;Fp,e^*R/,(0)F^,, 

All Ai2 
A21 A22 



p,ev 



Ht?H 



Fp,e«*R;,(0)*^F 



p,ev 



(14) 



where 77 = rt(l) = rt(— 1), since rt{m) = rt{—m) for 
the WSS Gaussian process h{n,T). If f] is pre-known — for 
example, when the maximum Doppler, denoted as fa, can be 
measured, according to the Jakes' model, (1) = Jo(27r/dTs), 
where Jo( ) is the zeroth order Bessel function of the first kind 
— we can ehminate the effect of Doppler from Rcon.p by 

All Ai2/r? 
A2i/?7 A22 



Rc 



(15) 



However, when rj can not be obtained in advance, we can 
approximately evaluate 77 by 



IIA21I 



2 ||2 

IIAiillp, -|- IIA22IIP' 



(16) 



Noting when the pilots are equispaced and P is a factor of 
N, which is an optimal case [12], (16) is accurate. In fact, at 
this circumstance, F^ev(od)^p.et^(od) = P'^l, therefore 

L-l 

||Aii||^ = ||A22||^ = P'^a* (17) 



||Ai2||^ = ||A2i||^ = v'P'T.''' 



(18) 



=0 



From (17)(18), (16) is an accurate estimator of rj. 

Denoting Fp = FJ^^, (Fp,e„*)^ 



Rc 



, (15) is rewritten into 
FR^(0)f'' (19) 



and with the noise component of (12), (19) is 

Rcon,p=Rcon,p + (^Ihp = FR^O)^'^ + a%P (20) 

Now (20) is of the standard form to apply ESPRIT algo- 
rithm, which we will briefly discuss in the following. Firstly, 

the eigenvalue decomposition (EVD) is applied to Rcon,p- 
Then L dominant eigenvalues are distinguished and the corre- 
sponding eigenvectors are collected in a matrix denoted as U. 
As span{lJ) = span{Fp), there exists a nonsingular matrix 
T such that U = FpT. Vertically splitting U into two parts, 
i.e., U = [Ujp,Uj^]^, it follows that U„j, = Fp,e«T and 
= Fp^ei;*T. Accordingly we have 

Vdw = U„j,T-i*T = U„pQ (21) 

where Q = T~^$T is similar with in other words, they 
have the common eigenvalues. We can solve Q from (21) by 
the linear least-squares (LS) criterion or the non-linear total 
least-squares (TLS) criterion for a better solution. When Q is 
obtained, its eigenvalues, denoted as A;, are calculated through 
EVD. Then the tap delays are estimated as 

arg{X;)N 



n 



2ttv 



,1 = 0, 



,L-1 



(22) 



where arg{X) denotes the phase angle of A in the range [0, 2n). 
Finally, as the phase angle wraps around with a period of 27r, 
the multipath delay is uniquely identified only when r < N/ v, 
where r is normalized by the sampling period T. 

C. Estimation of the Number of Significant Paths 

In practical applications, the auto-correlation matrix Rcon,p 
in (12) is unattainable. With a finite number of received 
symbols, the sample auto-correlation matrix is obtained as 



Rl, 



,p(n)Hf„„,p(n) (23) 



where Nt is the number of sample OFDMA symbols. Rcon,p 
is an asymptotic unbiased consistent estimator of Rcon.p- 

Since 'Rcon,p is affected by the Doppler, it is not used in 
the estimation of the number of significant paths. Instead, the 
auto-correlation matrix of Hp in (5) is adopted. With (14), it 
can be easily obtained through R, 



■con,p 



as 



Rp = (Aii+A22)/2 



(24) 



where we note *Rft(0)*^ = R/i(0). 

When R',o„ p is available, the sample auto-correlation ma- 
trix RJ^ is obtained through (24) accordingly, with which, then, 
MDL is performed to estimate the number of significant paths, 
denoted as L. 

D. Channel Interpolation 

After estimating the multipath delays, further modifications 
are made to enhance the performance. Due to many potential 
factors, such as the inaccurate estimation of rj, finite obser- 
vations of CFR and additive noise, ESPRIT is impaired: the 
estimated multipath delays fluctuate around their true values 



within a certain range, which results in an incomplete path 
subspace spanned. In addition, since the multipath delays are 
non-sample-spaced, all paths leak their power to the samples 
nearby. And the closer a sample gets to a path, the stronger it 
is influenced by the path. 

Considering the fluctuating estimated delays as well as the 
power leakage, we prefer to broaden the "observation window" 
around each estimated path delay to capture most of its power, 
i.e., for each estimated path f;, a capturing window, denoted 
as Wfi = {[ti\ — /5, + is included, where /3 is 

a predefined parameter such that 2/3 L < P. Through the 
simulation, we find /3 < 5 is sufficient for a wide range 
of SNR's and Doppler's. Consequently, the expanded set of 
estimated multipath delays If is 

If = f|Icp (25) 

where Tcp = {0, Lcp — 1} denotes the range of CP, hence 
limits the estimated multipath delays within it. 

The parameter /? requires carefully designing, since there 
exists a tradeoff: the larger f3 is, the more energy of paths is 
captured, meanwhile, the more noise is also introduced, and 
vice versa. 

The channel interpolators, denoted as Gev and God for even 
and odd symbols, respectively, are given as 

Gev{od) = Fd,ev(od)Fley(^od) (^6) 

where Fd ,,y(^od) is a matrix with the (fc, Z)-th element 
[Fd,ev{od)]k,i = e-^^'^'^*^'/^, where (j)k £ Id,ev{od) and f; e 
2f, respectively; similarly, Fp e„(o(;) is a matrix with the (fc, l)- 
th element Fp^et,(od) j,.; = e^J^Trekfi/Ar^ ^^^^.^ ^ Tp^ev(od) 
and fi e If, respectively; Id.ev(od) denotes the indexes of data 
tones allocated to the given user in the even (odd) OFDMA 
symbols. 

Finally, the channel estimation on the data tones is 

'tid,ev(od) — Gev{od)Flp,ev{od) (27) 

where flp_ev{od) denotes the LS channel estimation on the pilot 
tones of even (odd) OFDMA symbols. 

IV. Simulation Results 

The performance of the proposed channel estimation algo- 
ritiim is evaluated for a WiMAX system [1] with BW — 10 
MHz, N = 1024, Nused = 840, L^p = 128, /c = 3.5 
GHz. For simplicity, data tones are QPSK modulated, and no 
forward error coding is applied. ITU Vehicular A channels 
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Fig. L NMSE comparison for LL, PH and DP under different A^t's when 
fd = 200 Hz and Af^ch = 20. 

[13] is adopted, which consists of six individually faded taps 
with relative delays as [0, 310, 710, 1090, 1730, 2510] ns 
and average power as [0, -1, -9, -10, -15, -20] dB. The 
Jakes's spectrum [9] is applied to generate the Rayleigh fading 
channel. Besides, ideal synchronization is assumed, thus the 
delay of the first tap is always zero (tq = 0). 

The performance of the proposed pilot hopping based chan- 
nel estimation algorithm (PH) is compared with the local linear 
interpolation algorithm (LL) when uplink virtual MIMO is 
active. The algorithm in [5] which can be called the "doublet- 
pilots" algorithm (DP) is also evaluated as a benchmark, and 
for it, MIMO is inactive. For LL, the CFR of a certain data tone 
is simply equal to the arithmetic average of the two pilot tones 
within the same tile, since the irregular distributed pilot pattern 
prohibits the global interpolation. Since the second symbol in 
a tile contains no pilots, PH is adjusted to accommodate this 
situation: the second symbol in each tile is removed and CFR 
of this symbol is linearly interpolated from the first and third 
symbols within the same tile. There are three main factors 
influence the performance of PH and DP dominantly: the 
number of observed OFDMA symbols (iVj), Doppler {jd), and 
the number of allocated subchannels Q (Nsch)- 

First, different iV^'s are evaluated to examine the conver- 
gence of PH. In Fig.l, three values of Nt, i.e., 10ms, 20ms 
and 40ms, which are 96, 192 and 387 OFDMA symbols, cor- 
respondingly, are plotted when fd — 200 Hz and Nsch = 20. 
Obviously, the performance of LL is independent of Nt, and it 

'Eacii subcliannel consists of six tiles random distributed in frequency 
domain. 
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Fig. 2. NMSE comparison for LL, PH and DP under different /^j's wlien 
Aft = 192 and N^ch = 20. 
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levels off at high SNR regime. For PH and DP, no significant 
error floor can be observed. When Nt ^ 10 ms, PH at least has 
a 2.5dB gain over LL at low SNR regime, and far better than 
LL at high SNR regime. Moreover, PH is about 2dB worse 
than DP for all Nt's. However, it is worth noting that for the 
same Nt, the number of pilots available for DP is two times 
for PH, since MIMO is active for PH but not for DP. 

Fig. 2 plots the performances of LL, PH and DP under 
different Doppler's, when Nt — 192 and Ngch = 20. For LL, 
the performance gets worse along the increasing of f^, and 
levels off at high SNR regime for all f/s. However, for PH 
and DP, no significant error floor can be observed for lower 
/d's. Although PH also levels off when fd is high, it still 
outperforms LL over lOdB at high SNR regime. Moreover, 
when fd is lower, the performance difference between PH and 
DP is subtle, and when fd is higher, the NMSE difference is 
about 5dB at high SNR regime. 

Finally, different Nsch'& are evaluated in Fig. 3. From the 
figure, it is obvious that the performance difference between 
PH and DP is subtle when Nsch > 20, and both of them are far 
better than LL. When Nsch = 10, PH outperforms LL about 
2dB at low SNR regime and over 8dB at high SNR regime, 
meanwhile, it is worse than DP for about IdB at low SNR 
regime and 6dB at high SNR regime. 

From simulations, we find that PH outperforms LL over 
a wide range of SNR's and Doppler's when the number 
of observed OFDMA symbols and the number of allocated 
subchannels are not too small, i.e., Nt > 96 and Ngch > 10- 

V. Conclusion 

In this paper, we propose an ESPRIT-based channel esti- 
mation algorithm applicable to uplink OFDMA by exploiting 
pilot hopping. Through a very simple pilot hopping scheme, 
the shift-invariance property based on which ESPRIT is capa- 
ble is acquired. Since this special property is derived from a 
pair of contiguous OFDMA symbols, no special pilot pattern, 
e.g., the doublet pilots in [5], is indispensable within one 
OFDMA symbol. Hence, the proposed algorithm increases 



the spectrum efficiency and eases the design of the pilot 
pattern. Estimating the auto-correlation matrix over numbers 
of OFDMA symbols, the proposed algorithm outperforms 
the linear local interpolator within a wide range of SNR's 
and Doppler's. Besides, low rank adaptive filter [14] can be 
integrated to abate the estimation latency. 
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